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Abstract 

We consider the general situation of a compact relativistic body with anisotropic pres- 
sures in the presence of the electromagnetic field. The equation of state for the matter 
distribution is linear and may be applied to strange stars with quark matter. Three 
classes of new exact solutions are found to the Einstein-Maxwell system. This is 
achieved by specifying a particular form for one of the gravitational potentials and 
the electric field intensity. We can regain anisotropic and isotropic models from our 
general class of solution. A physical analysis indicates that the charged solutions de- 
scribe realistic compact spheres with anisotropic matter distribution. The equation of 
state is consistent with dark energy stars and charged quark matter distributions. The 
masses and central densities correspond to realistic stellar objects in the general case 
when anisotropy and charge are present. 
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1 Introduction 



Since the pioneering paper by Bowers and Liang [1] there have been extensive investiga- 
tions in the study of anisotropic relativistic matter distributions in general relativity to 
include the effects of spacetime curvature. The anisotropic interior spacetime matches 
to the Schwarzschild exterior model. The early work of Ruderman [2] showed that 
nuclear matter may be anisotropic in density ranges of 10 15 gem -3 where nuclear inter- 
actions need to be treated relativistically. Note that conventional celestial bodies are 
not composed purely of perfect fluids so that radial pressures are different from tan- 
gential pressures. Anisotropy can be introduced by the existence of a solid stellar core 
or by the presence of a type 3A superfluid as indicated by Kippenhahn and Weigert 
[3]. Different kinds of phase transitions (Sokolov [4]) or pion condensation (Sawyer 
[5]) can generate anisotropy. Binney and Tremaine [6] have considered anisotropies in 
spherical galaxies in the context of Newtonian gravitational theory. Herrera and Santos 
[7] studied the effects of slow rotation in stars and Letelier [8] analysed the mixture 
of two gases, such as ionized hydrogen and electrons, in a framework of a relativis- 
tic anisotropic fluid. Weber [9] showed that strong magnetic fields serve as a vehicle 
for generating anisotropic pressures inside a compact sphere. Some recent anisotropic 
models for compact self-gravitating objects with strange matter include the results of 
Lobo [10] and Sharma and Maharaj [11] with a barotropic equation of state. There- 
fore the study of anisotropic fluid spheres in static spherically symmetric spacetimes is 
important in relativistic astrophysics. 

In recent years there have been several investigations of the Einstein-Maxwell sys- 
tem of equations for static spherically symmetric gravitational fields usually with 
isotropic pressures to include the effects of the electromagnetic field. The interior 
spacetime must match at the boundary to the Reissner- Nordstrom exterior model. The 
models generated can be used to describe charged relativistic bodies in strong gravi- 
tational fields such as neutron stars. Many exact solutions have been given by Ivanov 
[12] and Thirukkanesh and Maharaj [13] which satisfy the conditions for a physically 
acceptable charged relativistic sphere. Charged spheroidal stars have been studied ex- 
tensively by Komathiraj and Maharaj [14], Sharma et al [15], Patel and Koppor [16], 
Tikekar and Singh [17] and Gupta and Kumar [18]. These charged spheroidal models 
contain uncharged neutron stars in the relevant limit and are consequently relevant in 
the description of dense astrophysical objects. We point out the particular detailed 
studies of Sharma et al [19] in cold compact objects, Sharma and Mukherjee [20] anal- 
ysis of strange matter and binary pulsars, and Sharma and Mukherjee [21] analysis 
of quark-diquark mixtures in equilibrium in the presence of the electromagnetic field. 
Charged relativistic matter is also relevant in modeling core-envelope stellar system 
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as shown in the treatments of Thomas et al [22], Tikekar and Thomas [23] and Paul 
and Tikekar [24] in which the stellar core is an isotropic fluid surrounded by a layer of 
anisotropic fluid. Consequently the study of charged fluid spheres in static spherically 
symmetric spacetimes is of significance in relativistic astrophysics. 

From the above motivation it is clear that both anisotropy and the electromag- 
netic field are important in astrophysical processes. However previous treatments have 
largely considered either anisotropy or electromagnetic field separately. The intention 
of this paper is to provide a general framework that admits the possibility of tangential 
pressures with a nonvanishing electric field intensity. We believe that this approach 
will allow for a richer family of solutions to the Einstein-Maxwell field equations and 
possibly provide a deeper insight into the behaviour of the gravitational field. On 
physical grounds we impose a barotropic equation of state which is linear, that relates 
the radial pressure to the energy density and allows for the existence of strange matter. 
Our general model will contain strange matter solutions found previously. In this re- 
gard we mention the following recent works on strange stars. Mak and Harko [25] and 
Komathiraj and Maharaj [26] found analytical models in the MIT bag model (Witten 
[27]) with a strange matter equation of state in the presence of an electromagnetic field. 
Sharma and Maharaj [11] generated a class of exact solutions which can be applied to 
strange stars with quark matter for neutral anisotropic matter. Lobo [10] found sta- 
ble dark energy stars which generalise the gravastar model governed by a dark energy 
equation of state. 

The objective of this treatment is to generate exact solutions to the Einstein- 
Maxwell system, with linear equation of state, that may be utilised to describe a 
charged anisotropic relativistic body. In Section 2, we express the Einstein-Maxwell 
system as a new system of differential equations using a coordinate transformation, 
and then write the system in another form which is easier to analyse. Three classes of 
new exact solutions to the Einstein-Maxwell system are found in Section 3 in terms of 
simple elementary functions. We show that particular uncharged anisotropic strange 
stars found in the past are contained in our general family of solutions. In Section 
4, we show that the solutions are physically admissible and plot the matter variables 
for particular parameter values. We generate values for the mass and central density 
in Section 5 for charged and uncharged matter. This analysis extends the treatment 
of Sharma and Maharaj [11] to include charge, and confirms that the exact solutions 
found are physically reasonable. Some concluding remarks are made in Section 6. 
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2 The field equations 



Our intention is to model the interior of a dense star. On physical grounds it is necessary 
for the gravitational field to be static and spherically symmetric. Consequently, we 
assume that the interior of a spherically symmetric star is described by the line element 

dg 2 = _ e 2u(r) dt 2 + e 2X(r) dr 2 + ^2 + ^2 ^2) ^ 

in Schwarzschild coordinates (x a ) = (t,r,6,<p). We take the energy momentum tensor 
for an anisotropic charged imperfect fluid sphere to be of the form 

Tij = diag(-p - ^E 2 ,p r - ^E 2 ,p t + l -E\p t + ^ 2 ), (2) 

where p is the energy density, p r is the radial pressure, p t is the tangential pressure and 
E is the electric field intensity. These quantities are measured relative to the comoving 
fluid velocitry u l = e~ u 5 % Q. For the line element ([T]) and matter distribution §2§ the 
Einstein field equations can be expressed as 

i[r(l-e-)]' = P+ \e\ (3) 

e- a ( v » + S' + --v'X 1 --) = P, + \e\ (5) 
\ r r J 2 

a = ^e-\r 2 E)>, (6) 

where primes denote differentiation with respect to r and a is the proper charge density. 
In the field equations ((31)-®, we are using units where the coupling constant = 1 
and the speed of light c = 1. The system of equations ((31)-® governs the behaviour of 
the gravitational field for an anisotropic charged imperfect fluid. Note that the system 
®-® becomes 

^Hl-O]' = P, (7) 

-\(l-e- 2X ) + — e- 2 " = p, (8) 

e- 2 *(v" + v> 2 + --v>\'--) = p, (9) 
\ r r J 

for matter distributions with isotropic pressures (p r = pt) in the absence of charge 
(E = 0). 

The mass contained within a radius r of the sphere is defined as 
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r 



m ( r ) = 2 / u;2 p( u; )^ u; - (10) 



A different, but equivalent, form of the field equations is obtained if we introduce a 
new independent variable x, and define functions y and Z, as follows 



x = Cr 2 , Z(x) = e~ 2X{r) and A 2 y 2 



which was first suggested by Durgapal and Bannerji [28]. Then the line element ([I]) 
becomes 

ds 2 = -A 2 y 2 dt 2 + — ^dx 2 + (d6 2 + sin 2 Qd<\> 2 ) . (12) 
ACxZ C 

In (fTTj) and (fl2l . the quantities A and C are arbitrary constants. Under the transfor- 
mation (fTTj) . the system ©-© becomes 

l-Z 



x 



2Z 
- 1 



p E 1 
C + 2C 1 



■ y 



y x 
V , 



AxZ- + (AZ + 2xZ) 



p l _& 

C 2C" 
^ E 2 
VZ ~ C + 2C> 

<7 = ~( XE + E 



(13) 
(14) 

(15) 

(16) 



where dots denote differentiation with respect to the variable x. The mass function 
(TTOT) becomes 



m[x) 



4C73/2 



\fwp{w)d 



w. 



(17) 



in terms of the new variables in (fTTj) . 

For a physically realistic relativistic star we expect that the matter distribution 
should satisfy a barotropic equation of state p r = p r (p)- For our purposes we assume 
the linear equation of state 

p r = ap — P, (18) 

where a and P are constants. Then it is possible to write the system ( 1T3T) - ( TT6T) in the 
simpler form 

1 - Z 



P_ 

C 

Pr 
Pt 

A 



£2 
2C 

C 



x 

ap — P, 

Pr + A, 

ACxZ^- + 2C 

y 

(1 + ha 



■ E 2 
2Z , 

2C" 



xZ + 



AZ 



1 + a 



y 
y 



(1 + a) 

1 - z 



cz-^^l + 



2P 



1 + a) 



1 



x 



(1 + a) 



2aZ + AZ^ + ^ 
y c . 



A—(xE + E)' 

x 



(19) 
(20) 
(21) 



(22) 
(23) 
(24) 



where the quantity A = pt~p r is the measure of anisotropy in this model. In the system 
CEHD-dMD, there are eight independent variables (p,p r ,p t , A, E,cr,y, Z) and only six 
independent equations. This suggests that it is possible to specify two of the quantities 
involved in the integration process. The resultant system will remain highly nonlinear 
but it may be possible to generate exact solutions. 



3 Generating exact models 

We must make physically reasonable choices for any two of the independent variables 
and then solve the system (|T9"l) - (j24"j) to generate exact models. In this paper, we choose 
forms for the gravitational potential Z and electric field intensity E. We make the 
specific choices 

1 + (a - b)x 

z = — 7- , 25 

1 + ax 

_ k(3 + ax) 

C ~ (1 + ax) 2 ' 1 j 

where a, b and k are real constants. The gravitational potential Z is regular at the origin 
and well behaved in the stellar interior for a wide range of values for the parameters 
a and b. The electric field intensity is continuous, bounded and a decreasing function 
from the origin to the boundary of the sphere. Therefore the forms chosen in fT25l) - (!26l) 
are physically reasonable. On substituting ff25]) and fl26l) in ff23l) we obtain 

y (1 + a)b ab 



y 4 [1 + (a - b)x] 2(1 + ax) [1 + (a - b)x] 
(3(1 + ax) (1 + a)fc(3 + ax) 



AC [1 + (a - b)x] 8(1 + ax) [1 + (a - b)x) ' 

which is a linear equation in the gravitational potential y. For the integration of 
equation ff27j) it is convenient to consider three cases: b = 0, a = b and a ^ b. 



3.1 The case 6 = 



When 6 = 0, (1271) becomes 



y 
y 



AC 



with solution 



-k(l + a) 

D(l + ax) a exp 



f a)k(3 + ax) 
8(1 + ax) 2 

2fc(l + a) 



a(l + ax) 



(3x 
AC 



(28) 



(29) 



where D is the constant of integration. We observe that p = — ^- for this case which 
we do not consider further to avoid negative energy densities. 
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3.2 The case a = b 

When a = b, (J2"7I) becomes 

y (1 + a)a aa + ax) (1 + a)/c(3 + ax) . . 

y = 4 + 2(1 + ax) 46 8(1 + ax) ' ^ ^ 

On integrating (130]) we get 

2aa—(l + ct)k 

y = D{l + ax) a exp[F(x)], (31) 



where 



F(x) = [-K7(l + a) - 2(3 + a(26(l + a) - 
oG 



and D is the constant of integration. Then we can generate an exact model for the 
system (fl9>(|2lj) as follows 

e 2A = 1 + ax, (32) 

e 2u = A 2 D 2( 1 + ax ) 2aa -^ 1+a) exp[2F(g)], (33) 
_P (2a - fc) (3 + ax) 

C 2 (1 + ax) 2 ' 1 j 

p r = ap — ft, (35) 

p t = p r + A, (36) 

A = 1 (6 2 [fc 2 (l + a) 2 x(3 + ax) 2 

166(1 + ax) d 1 L 

+4a 2 x(3 - 8a + 9a 2 + a 2 (l + a) 2 x 2 + 2ax(2 + 3a + 3a 2 )) 

-4Jfe(12 + a 3 (l + a)V + a 2 x 2 (7 + 9a + 6a 2 ) + ax(12 + 5a + 9a 2 ))] 

-4Cx(l + ax) 2 [(l + a)(2a 2 x - 3k) - aft(k(l + a) - 6a - 4)] 

+4/3 2 x(l + ax) 4 } , (37) 
E 2 _ k(3 + ax) 

C ~ (1 + ax) 2 ' m 
in terms of elementary functions. 

The solution f[3"2"j) - f[3"5j) may be used to model a charged anisotropic star with a 
linear equation of state. In this case the mass function is 

(2a-A;)x 3 / 2 . , 

which is similar to forms used in other investigations. The gravitational potentials and 
matter variables are continuous and well behaved in the stellar interior. Note that 
when k = the model fl32|) - fl38|) reduces to a solution for uncharged anisotropic stars. 
Equation ([3"T]) yields 



A = — — -r {C 2 a 2 x [3 - 8a + 9a 2 + a 2 (l + a) 2 x 2 + 2ax(2 + 3a + 3a 

4o ( 1 ~ h~ (XX j 

-2Cx(l + ax) 2 [(l + a)a 2 x + a(3(3a + 2)] + (3 2 x(l + ax) 4 } (40) 
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when k = so that the model is necessarily anisotropic with A 7^ in general even in 
the simpler case of uncharged matter. Some treatments of the physical properties of 
anisotropic spheres in general relativity include the investigations of Dev and Gleiser 
[29, 30], Mak and Harko [31, 32], Chaisi and Maharaj [33, 34] and Maharaj and Chaisi 
[35] with A ^ 0. 



3.3 The case a ^ b 

On integrating (j2"7j) we get 



y = D(l + ax) m [l + (a - b)xf exp 



-a(3x 



_AC(a-b)_ 

where D is the constant of integration, and m and n are given by 
2a6 - (1 + a)k 



(41) 



m 



n 



4b 



[2a 2 C{k{l + a) - 2ab) - abC{5k(l + a) - 26(1 + 5a)) 



8bC(a - b) 2 
+b 2 {3kC(l + a) - 26(7(1 + 3a) + 2(3)} . 

Then we can generate an exact model for the system (fT9"]) - ([24"]) in the form 
1 + ax 



e 2X = . 

1 + (a — b)x' 

e 2u = A 2 D 2 {l + ax) 2m [l + {a-b)x] 2n exp 

p (26 -k) {3 + ax) 
C 2 (1 + ax) 2 ' 

p r = ap — j3, 



-a/3x 



2C(a - b) 



p t = p r + A, 
A ~ bC 



(42) 
(43) 

(44) 

(45) 
(46) 



6(7(1 + 5a) 



+ 



2(3 



(1 + ax) (l + a)(l + ax) 2 (1 + a) 



Cx[l + (a - b)x] 
(1 + ax) 



, a 2 m(m — 1) 
4 + 

(1 + ax) 2 



2a(a — b)mn 



1 + ax)[l + (a — b)x] 



(a — b) 2 n(n — 1) \ 2a(3[a{m + n) [1 + (a — b)x) — bn) 



a 2 (3 2 



[l + (a-6)x] 2 J (a-b)C(l + ax)[l + (a-b)x] AC 2 (a-b) 2 
4[1 + ax{2 + (a - b)x)\ - 6(5 + a)x 



E 2 
C 



2(a - 6)(1 + a)(l + ax) 3 [l + (a - b)x] 
[-Ab 2 Cn + a 3 a;(-4C(m + n) + (3x) + a 2 (4C(m + n){2bx - 1) + 0(2 - 6x)x) 

+a(-46 2 C(m + n)x + (3 + b(ACm + 8Cn - (3x))] , (47) 
k(3 + ax) 



[1 + ax) 2 



(4? 
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in terms of elementary functions. 

Therefore we have generated a second class of solutions (|42l) - (|48|) that models a 
charged anisotropic star with a linear equation of state. The mass function has the 
form 

The form of the mass function (I49p represents an energy density which is monoton- 
ically decreasing in the stellar interior and remains finite at the centre x = 0. It is 
physically reasonable and has been used in the past to study the properties of isotropic 
fluid spheres: Matese and Whitman [36] generated equilibrium configurations in gen- 
eral relativity, Finch and Skea [37] studied neutron star models and Mak and Harko 
[32] analysed anisotropic relativistic stars with this form of mass function. Lobo [10] 
demonstrated that (|49|) is consistent with stable dark energy stars which generalises 
the gravastar model of Mazur and Mottola [38]. It was then shown that large stability 
regions exist close to the event horizon thereby making it difficult to distinguish dark 
energy stars from black holes. Sharma and Maharaj [11] found a new class of exact 
solutions to Einstein equations that can be applied to strange stars with quark matter 
with this mass distribution. Consequently the mass function (j4"9|) is of astrophysical 
importance in the description of compact objects. 

It is interesting to observe that for particular parameter values we can regain un- 
charged anisotropic and isotropic models (k = 0) from our general solution (|42p - (|48p . 
We regain the following particular cases of physical interest: 

3.3.1 Sharma and Maharaj model 

If we set (3 = ap s then 

p r = a(p - p s ), 

where p s is the density at the surface r = s. Thus we regain the equation of state of 
Sharma and Maharaj [11]. Then by setting C = 1 and A 2 D 2 = B we find that the line 
element is of the form 

ds 2 = -B(l + ar 2 ) a [l + (a - b)r 2 \< exp ( ) dt 2 

\Z{a — b) J 

i + ax^ 2 2 s . n2 2 
1 + (a - b)x 

where 

5aba — 2a 2 a — 3b 2 a + ab — b 2 + b(3 
7= 2{a~b) 2 " 

The line element (1501) corresponds to the uncharged anisotropic model of Sharma and 
Maharaj [11]. They showed that this solution may be used to describe compact objects 
such as strange stars with a linear equation of state with quark matter. 
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3.3.2 Lobo model 



If we set (3 = then 

p r = ap 

and we regain the equation of state studied by Lobo [10]. Then on setting a = 2b, C = 1 
and A 2 D 2 = 1 we generate the metric 

ds 2 = -(1 + br 2 )^' 2 {l + 2br 2 ) a dt 2 + ( \ + + 2 ^ dr 2 + r 2 (d6 2 + sin 2 Odtf). (51) 

The line element (l5Tj) corresponds to the uncharged anisotropic model of Lobo [10]. 
We point out that the line element ( 1ST]) serves as an interior solution with a < — | 
which may be matched to an exterior Schwarzschild solution in a model for dark energy 
stars. Lobo [10] proved that stability regions exist for dark energy stars by selecting 
particular values of a in a graphical analysis. 

3.3.3 Isotropic models 

In general A ^ so that the model remains anisotropic. However for particular 
parameter values we can show that A = in the relevant limit in the general solution 
(S2D-(@HD • If we set a = and 6=1 then we obtain 

a 

m = - 

n = -Lp-(l + 3a)C] 

A = 4C(1- x) [P ~ 3(1 + a)C] [(3 - {1 + 3a)C] - (52) 

Two different cases arise as a consequence of (1521) if we set A = 0. 

In the first case we observe that when (3 = and a = — 1 then A = 0. The equation 
of state becomes p r (= Pt) = —p- In this case the line element becomes 

ds 2 = - (l - dt 2 + (l - dr 2 + r 2 {dB 2 + sin 2 Odcf) 2 ), (53) 

where we have set A = D = 1 and C = . The metric (153]) corresponds to the familiar 
isotropic uncharged de Sitter model. 

In the second case we see that when (3 = and a = — | then A = 0. The equation 
of state becomes p r {— Pt) = ~\p- 111 this case the line element becomes 

E 2 



ds 2 = -A 2 dt 2 + ( 1 - 1- ) dr 2 + r 2 (d# 2 + sin 2 ^0 2 ), (54) 



where we have set D = 1 and C = The metric (|54p corresponds to the well known 
isotropic uncharged Einstein model. 
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4 Physical analysis 



The solutions found in this paper may be connected to the Einstein-Maxwell equations 
for the exterior of our source. We need to match the Reissner-Nordstrom exterior 
spacetime 

ds 2 = - ( 1 + -V eft 2 + 1 + -^r dr 2 + r 2 (d# 2 + sin 2 ddcj) 2 ) 

to the interior spacetime (TTT) across the boundary r = R. This generates the conditions 



2 



2M Q 
~R + R? 



A 2 y 2 (CR 2 ) 



2M Q 2 \ 1 1 + aCR 2 



R R? J l + (a-b)CR 2 

which relates the constants a, b,A,C,D,at and (3. This demonstrates that the continuity 
of the metric coefficients across the boundary of the star r = R is easily satisfied as 
there are sufficient number of free parameters. If there is a surface layer of charge then 
the pressure may be nonzero which would place restrictions on the function v through 
the matching conditions at the boundary. However the number of free parameters 
available easily satisfies the necessary conditions that arise for a particular model under 
investigation. 

We now briefly consider the physical behaviour of the models generated in Section 
3 for the case a ^ b. From the explicit forms (1421) and (1431) we can easily see that 
the gravitational potentials e 2v and e 2A are continuous, well behaved and nonsingular 
at the origin. The energy density p is continuous and monotonically decreasing from 
the centre to the boundary of the star, which is a necessary condition for a realistic 
model. The radial pressure p r also has the same feature because p and p r are linked 
by a linear equation of state. The tangential pressure p t is also nonsingular at the 
origin and continuous for a wide range of the parameters a, b and k. To maintain the 
usual casuality condition we must place the restriction that < a < 1 if we require 
-J^ < 1. However note that our models do allow for a < in the case of anisotropic 
dark energy stars. The form chosen for electric field intensity E is physically reasonable 
and describes a decreasing function . 

With the help of a particular example we can demonstrate the above features graph- 
ically. Figures 1-4 represent the energy density, the radial pressure, the tangential 
pressure and the electric field intensity, respectively. To plot the graphs we choose the 
parameters a = 3, b = 2.15, a = 0.33,/? = ap s = 0.198, C = 1 and k = 0.2, where p s is 
the density at the boundary r = s = 1.157. Note that our choice of a = 0.33 ensures 
that both the radial pressure and the tangential pressure for the neutral sphere vanish 
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at the boundary. We observe from Figures 1-4 that the matter variables p } p r ,Pt and 
E have the appropriate features to describe a compact relativistic sphere. Solid lines 
represent uncharged matter and dashed lines include the effect of charge in Figures 
1-3. We observe that the effect of E is to produce lower values for p,p r and p t when 
compared to the case of uncharged matter. In Figure 5 we have plotted the measure 
of anisotropy A for the same parameter values used above. Note that the effect of the 
electromagnetic field is to increase the magnitude of A which affects the behaviour of 
Pt- 





1.0 1.2 



Figure 1: Energy density. 





0.2 0.4 0.6 0.8 1.0 1.2 



Figure 2: Radial pressure. 
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Figure 3: Tangential pressure. 



Figure 4: Electric field intensity. 

5 Stellar structure 

In this section we show that the solutions generated in this paper can be used to 
describe realistic compact objects. In particular we seek to compare our results with 
those of Sharma and Maharaj [11] since they regain values for the stellar mass agreeable 
with observations. To achieve consistency with Sharma and Maharaj [11] we introduce 
the transformations 

a = aR 2 , b = bR 2 , (3 = (3R 2 , k = kR 2 . 

Under these transformations the energy density becomes 

(2b — k)(3 + ay) 



P 



2R 2 (1 + ay) 2 



(55) 



O.d 



0.8 _1.Q_ _ 1.2 



-0.1 : 
-0.2 1 
-0.3 ; 
-0.4 ; 
-0.5 ; 
-0.6 - 



E*0 



Figure 5: Measure of anisotropy. 
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and the mass contained within a radius s has the form 

(2b-~ky/R* 
= 4{i + ~asy&y (56) 

where we have set C — 1 and y = When k = (or £7 = 0), (I55I) and (|56|) reduce 
to the expressions of Sharma and Maharaj [11]: 

5(3 + ay) M _ fe 3 /^ 2 , - 

P i? 2 (l + m/) 2 ' 2(1 + as 2 /R 2 ) 1 J 

which gives the density p and mass M of an uncharged star of radius s. 

If we choose a = 53.34, b = 54.34, R = 43.245 km and s = 7.07 km then we 
can produce an uncharged model (k = 0) with mass M = 1.433Mq and central 
density p c = 4.672 x 10 15 gcm -3 . The corresponding value of a = 0.437 is obtained 
by requiring that the anisotropy vanishes at the boundary. To simplify comparison 
with Sharma and Maharaj [11] we have used the same values of a, b, R and s; however 
our value for a is a correction. It should be noted that these results are consistent 
with the equation of state for strange matter formulated by Dey et al [39] . This has 
astrophysical significance as their model has been used to describe the X-ray binary 
pulsar SAX J1808. 4-3658. When the charge is nonzero we set k = 37.403 and then 
we obtain the mass M = 0.940Mq and central density p c = 3.064 x 10 15 gcm -3 . The 
values for M and p c generalise the figures of Sharma and Maharaj [11] to include the 
effect of the electromagnetic field. Choosing different set of values for the parameters 
will produce different results as shown in Table 1. Note that the values presented in 
Table 1 correspond to a star of radius s = 7.07km. The value of k = 37.403 is selected, 
in generating Table 1, so that the density and mass of the Sharma and Maharaj [11] 
analysis is regained for uncharged matter. Furthermore, the value of k = 37.403 with 
E — generates a star of mass 1.433Mq which is the same as the strange star model 
of Dey et al [39]. With this value of k we find that the star has mass 0.940Mq in 
the presence of charge so that the stellar core has a lower density which represents a 
weaker field. This is consistent as the effect of the electromagnetic field is repulsive. 

We observe that the values for the mass in the presence of charge (E ^ 0) is always 
less than the uncharged case. The central density of the charged sphere is also less 
than the uncharged case. Sharma and Maharaj [11] showed that anisotropy affects the 
mass and central densities of massive objects. We have shown that the inclusion of the 
electromagnetic field also affects M and p c . Both anisotropy and charge are physical 
quantities that affect the range of degenerate states in our model. For the calculation 
of mass and central density we have set s = 7.07 km, R = 43.245 km, k = 37.403 and 
p s = 1.17119 x 10 15 gcm -3 for the uncharged case. 
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b 


a 


a 


Pc 


M 


Pc 


M 








(x 10 15 gcm~ 3 ) 


(Mq) 


(xlO 15 gem -3 ) 


(Mq) 








E = 


E = 


E ^ 


E^O 


30 


23.681 


0.401 


2.579 


1.175 


0.971 


0.443 


40 


36.346 


0.400 


3.439 


1.298 


1.831 


0.691 


50 


48.307 


0.424 


4.298 


1.396 


2.691 


0.874 


54.34 


53.340 


0.437 


4.671 


1.433 


3.064 


0.940 


60 


59.788 


0.457 


5.158 


1.477 


3.550 


1.017 


70 


70.920 


0.495 


6.017 


1.546 


4.410 


1.133 


80 


81.786 


0.537 


6.877 


1.606 


5.269 


1.231 


90 


92.442 


0.581 


7.737 


1.659 


6.129 


1.314 


100 


102.929 


0.627 


8.596 


1.705 


6.989 


1.386 


183 


186.163 


1.083 


15.730 


1.959 


14.124 


1.759 



Table 1: Central density and mass for different anisotropic stellar models for neutral 
and charged bodies 

6 Conclusion 

We have found a general framework for the Einstein-Maxwell system of equations with 
a linear equation of state for anisotropic matter distributions in the presence of the 
electromagnetic field. Three new classes of exact solutions have been generated to this 
system of nonlinear equations. We have shown that these classes of solutions satisfy the 
necessary physical requirements in the description of a charged compact objects with 
anisotropic matter distribution. We have demonstrated that our models yield stellar 
structures with masses and densities consistent with the Dey et al [38] and Sharma and 
Maharaj [11] models in the limit of vanishing charge. Therefore it is likely that our 
solutions may be helpful in the gravitational description of stellar bodies such as SAX 
J1808. 4-3658. The solutions obtained may be useful to model the interior of charged 
relativistic quark stars with anisotropic matter distribution. Our models contain the 
uncharged anisotropy models of Sharma and Maharaj [11] and Lobo [10] which de- 
scribe quark stars and strange matter stars. We believe that the general class of exact 
solutions found in this paper may assist in more detailed studies of relativistic stellar 
bodies and allows for different matter distributions because of the form of the linear 
equation of state chosen. 
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